arX1v:1904.02187v2 [hep-th] 30 Apr 2019 


Creating a ‘Traversable Wormhole 


Gary T. Horowitz’, Don Marolf’, Jorge E. Santos®, Diandian Wang’ 
* Department of Physics, University of California, Santa Barbara, CA 93106 


° Department of Applied Mathematics and Theoretical Physics, 
University of Cambridge, Wilberforce Road, Cambridge CB3 OWA, UK 


We argue that one can nucleate a traversable wormhole via a nonperturbative process in quan- 
tum gravity. To support this, we construct spacetimes in which there are instantons giving 
a finite probability for a test cosmic string to break and produce two particles on its ends. 
One should be able to replace the particles with small black holes with only small changes to 
the spacetime away from the horizons. The black holes are then created with their horizons 
identified, so this is an example of nucleating a wormhole. Unlike previous examples where the 
created black holes accelerate apart, in our case they remain essentially at rest. This is impor- 
tant since wormholes become harder and harder to make traversable as their mouths become 
widely separated, and since traversability can be destroyed by Unruh radiation. In our case, 
back-reaction from quantum fields can make the wormhole traversable. 
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1 Introduction 


Classically, there are no traversable wormholes. This is a consequence of topological censorship 
[I], which says that if the null energy condition is satisfied, any causal curve that starts and ends 
at infinity can be continuously deformed to a causal curve that stays in the asymptotic region. 
However, it has recently been shown that quantum matter fields can provide enough negative 
energy to allow some wormholes to become traversable. This was first shown for asymptotically 
anti-de Sitter wormholes using ideas from gauge/gravity duality [4| and later generalized 
to asymptotically flat wormholes [5 af 

While the wormhole may only be traversable for a limited time, in all of the examples that 
have been discussed so far the wormhole itself is eternal, i.e., it exists on all complete spacelike 
surfaces. This is not surprising since the quantum fields produce only small effects and there 
is a theorem which says that the topology of space cannot change in classical general relativity 
[8]. 

The topology of space can change through a nonperturbative quantum tunneling event. This 
tunneling process can be investigated in the path integral formulation of quantum gravity if a 
sum over all topologies is assumed in the integral. Instantons, classical solutions to Euclidean 
field equations, have been widely used to study non-perturbative quantum processes of this type 
[9]. For example, there is an instanton describing the pair creation of oppositely charged black 


'In these examples, it takes longer to go through the wormhole than to go around. As these papers explain, 
one expects this will be true for all traversable wormholes since it is required by the generalized second law 
and boundary causality (when there is a field theory dual). 


holes in a background electric or magnetic field [10]. This is directly analogous to the Schwinger 
pair creation of oppositely charged particles in an electric field. The black holes are created 
with their horizons identified, so this process creates a wormhole. In the small black hole limit, 
the instanton reduces to the standard one-particle description of the Schwinger effect as arising 
from a Euclidean solution where the black hole particle moves on a circle about the origin ; 
Several other analogous examples have also been found and studied [18]. 

If the mouths of the wormhole followed static trajectories, the path integral that computes 
quantum corrections to the wormhole would be of the form studied in [6]. So, as argued 
there, in a semiclassical approximation with the metric determined by the expectation value 
of the quantum stress tensor, and with appropriate boundary conditions, back-reaction from 
quantum fields would generically transform the solution into a traversable wormhole. Justifying 
this approximation typically requires that many degrees of freedom propagate in the throat [5]. 

Unfortunately, although the above instantons create black holes at rest, the background 
fields cause the black holes to quickly accelerate away from one another so that the analy- 
sis of [6] does not apply. Indeed, the black holes are separated by an acceleration horizon. 
Mathematically, this allows a Killing symmetry that generates both horizons and requires any 
quantum stress tensor to vanish when contracted with the horizon generators. This means that 
back-reaction from quantum fields will not make the wormhole traversable. Physically, one 
finds that the wormhole is created in an excited state dictated by the Unruh temperature de- 
fined by the acceleration horizons. So while the wormhole ground state may be traversable, the 
probability of finding such a state in the associated thermal ensemble is small due to both the 
tiny gap separating the energies of traversable and collapsing wormholes and the small number 
of traversable wormhole states relative to non-traversable black holes [4] [5]. If one tries to tune 
parameters to reduce the acceleration, the black holes are created farther and farther apart. 
Not only might one like to create wormhole mouths at finite locations, but a wide separation 
between the mouths makes the wormhole even more fragile so that the same problem remains. 

In particular, the asymptotically flat four-dimensional traversable wormholes of with 
radius ro and separation d have a large parameter q ~ ro/pianck that allows the gap Eyap ~ g/d 
between traversable and collapsing wormholes to be parametrically large compared with the 
temperature T = a set by the acceleration horizons. But the number of traversable wormhole 
states is small and the black hole entropy Spy = A/4¢ 


iani = TO is much larger than 


planck 
Egap/T ~ To/lplanck, sO non-traversable wormholes dominate the ensemble at all accelerations 
a ~ 1/d. On the other hand, in the appropriate background such traversable wormholes will 


dominate if we can find instantons with 
ad & Lplanck /T0- (1.1) 


We will indeed find backgrounds that admit instantons with vanishingly small accelerations, 
though ours will be asymptotically anti-de Sitter. 


The goal of this paper is to argue that one can produce traversable wormholes through a 
quantum tunneling event. As mentioned above, to achieve this we need to create the wormhole 
mouths a finite distance apart with arbitrarily small acceleration. We can think of this as 
creating a pair of black holes whose horizons are initially identified, with the understanding 
that back-reaction from quantum fields will then render the wormhole traversable. So long as 
the saddle-point black holes are near extremality and the theory contains massless fields that 
can propagate in the wormhole throat with appropriate boundary conditions, the wormhole 
will remain open for a long time M, [5] [6]. Fields of this type include massless bulk fields 
that allow s-waves (so unfortunately this fails for photons or gravitons), massless modes of 
fermions bound to magnetic field lines, and massless excitations on cosmic strings threading 
the wormhole. Under the proper boundary conditions, fields that propagate less well in the 
throat will still make the wormhole traversable but, if only such fields are present, the time the 
wormhole remains open is generally small and vanishes at extremality. 

For simplicity, we will model the black holes as point particles and consider a process in 
which a cosmic string breaks and nucleates the two particles on its ends. We will also neglect the 
backreaction of the cosmic string and created particlesP We will explicitly construct spacetimes 
in which a test string can nucleate two test particles with arbitrarily small acceleration. In 
effect, we find situations in which the standard circular Schwinger instanton mentioned above 
is distorted into a highly eccentric configuration whose extent in the Euclidean time direction 
is much larger than its width in space. One should be able to replace our particles with small 
black holes with negligible changes to the geometry elsewhere, and one may choose to identify 
their initial horizons as desired. If the black holes have the correct temperature to match the 
periodicity of the particle orbit in Euclidean time, then the full instanton will be smooth. But 
we may also consider instantons with conical singularities at the horizons as in [19]. 

In Section 2, we explain a subtlety with the instantons we use. In section 3 we derive the 
conditions on the spacetime in order to create nearly-static particles. The calculation of the 
action for our instantons is given in section 4. Section 5 provides a simple example of a star in 
AdS which satisfies our conditions, while section 6 describes a class of vacuum examples. The 
latter are obtained by deforming the metric on the boundary of global AdS. Section 7 concludes 
the paper with some final remarks. 


2 An Instanton Subtlety 


Before proceeding to the construction of suitable spacetimes, we explain a subtlety with the 
instanton we need for our problem. The energy required for the nucleation of black holes (or 
particles) is provided by the reduction in the potential energy of the system. For example, 
with a magnetic field background, the appearance of a magnetic monopole pair reduces the 


?For discussions of a gravitating cosmic string nucleating two black holes which accelerate apart, see [L7] [18]. 


electromagnetic potential energy; for the breaking of a cosmic string, the removal of a piece 
of string reduces the elastic potential energy. Unfortunately, the same feature that provides 
the required energy also provides a force for the created particles to accelerate away from 
each other. In order to create nearly static particles, we must somehow counterbalance this 
undesired acceleration. In this paper, we show that the curvature of spacetime can provide this 
counterbalance. 

One cannot create completely static particles since if the Lorentzian solution for particles 
attached to cosmic strings describes two static worldlines, the Euclidean solution will also 
consist of two worldlines which are independent of Euclidean time. So there will not be a slice 
of the Euclidean solution describing the initial unbroken cosmic string. Instead we want a 
family of instantons which describe the creation of particles with arbitrarily small acceleration, 
and thus finite but arbitrarily small Unruh temperature. 

It is instructive to think in terms of effective potentials. If the static Lorentzian solution 
has exactly the same energy as the unbroken cosmic string, then the effective potential has two 
degenerate minima (see Fig. [La). One minimum represents the unbroken string and the other 
represents the broken string with static particles on its ends. There is no “decay” from one 
minimum to the other in this double well potential. However, we will see that by adjusting 
the ratio of the string tension to the mass of the particles, one can lower the second minimum 
continuously. Now there are instantons which describe the nucleation of two particles on the 
ends of the string (see Fig. [1b). Assuming the original static solution is stable to small fluc- 
tuations (which will be the case in our examples) these particles will oscillate around the true 
minimum. We can now tune this acceleration to be as small as desired by adjusting the second 
minimum. The action for the instanton relative to the unbroken cosmic string, AS, remains 
finite in the limit that the minima become degenerate, but the decay probability goes to zero. 
This is because the the decay rate is given by 


T = |K]e7?^sS (2.1) 


where K is the contribution from small fluctuations about the instanton. This vanishes in the 
limit since the negative mode goes away. 

In this paper, we will calculate AS in the limit that the minima become degenerate and 
interpret it as providing the leading semiclassical decay rate when the minima are slightly 
separate. We will also refer to the static separation as the “separation between the created 
particles” with the understanding that this is the limiting case. As discussed in the introduction, 
our instantons will create traversable wormholes so long as the acceleration is sufficiently small, 
proportional to a power of Lplanck (See for the case of the wormholes of BF] We expect 


3 Alternatively, one can work with the strict double well potential. Since a state localized in the first minimum 
is a linear superposition of the symmetric and antisymmetric energy eigenstates, if we start with the unbroken 
cosmic string and wait, it will turn into a broken string with static particles on its end. At that time we can 
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(a) Degenerate minima (b) Nondegenerate minima 


Figure 1: (a) The red dot (unbroken string configuration) and black dot (broken string with 
static particles on the ends) are degenerate. In this case, the wavefunction of the true ground 
state (blue dashed curve) is a symmetrized state centered around both minima. (b) The red 
dot tunnels to the purple dot (broken string with particles momentarily at rest on ends), which 
then oscillates around the black dot (the true vacuum). As the static pair creation limit is 
approached, the black dot moves up towards the purple dot so the difference in the potential 
energy between the red dot and the black dot diminishes. 


K for such cases to be proportional to some positive power of the planck, and AS to diverge 


like an inverse power of Lplanck, SO the leading contribution to I as Cptanck — 0 is indeed a. 


3 Conditions on the Spacetime 


We are interested in static metrics that admit a two-dimensional totally geodesic sub-manifold. 
Totally geodesic submanifolds have vanishing extrinsic curvature and in particular are extremal 
surfaces. We may thus take this surface to be the worldsheet of our cosmic string in the case 
where it does not decay. So long as no external forces act transversely to this submanifold, for 
appropriate initial conditions the totally geodesic condition further allows us to compute the 
motion of particles using only the induced two-dimensional metric 
; 2 dr? 
dsźŹp f(r)dt* + aay (3.1) 
In particular, this will suffice to study the motion of particles created by the decay of our cosmic 
string. Furthermore, in the cases we will study below, the spacetime metric has an additional 
Zə symmetry r > =r. 
We would like to understand under what conditions one can find a nearly static particle 
pair created by cosmic string breaking in such backgrounds. We thus consider the following 


place an obstacle between the two masses to stop it from transitioning back into an unbroken cosmic string. 
This is another way to create a traversable wormhole. 


(a) Subsequent acceleration (b) Subsequent oscillation 


Figure 2: The process of string breaking. The green regions in the vertical and horizontal 
planes are the string worldsheet W in the Lorentzian and Euclidean sections respectively. The 
blue curve represents the Euclidean trajectory of the particle, while the red curves represent 
the worldlines of the physical particles after being created. (a) The traditional case where the 
two particles subsequently fly apart, being pulled by the broken pieces of string. (b) Our case 
where the particles are created close to local minima of an effective potential V(r), so that 
the subsequent motion is oscillatory with small acceleration. Note that the Euclidean solution 
in this case is not a symmetric circle but is elongated so the acceleration at the transition is 
smaller. 


action 


dz? dz? 
= — A= ; 
r=- ff dA—m m Gab Ta dà, (3.2) 


where the subscript L denotes Lorentzian signature, W is the worldsheet of the unbroken string 
or the union of both string worldsheets that arise from breaking the original cosmic string, dA 
is the proper area element on the worldsheet, A is a parameter along the boundary OW of W, 
u is the string tension, and m is the particle mass. In our set up, W has either no (inner) 
boundary or two boundaries, corresponding to an infinitely long string that is either unbroken 
or broken (see Fig. [2). 

We are interested in situations when the string is broken while preserving the Zə symme- 
try. It is convenient to subtract the action of the unbroken cosmic string and work with the 


as,-—» [™ f" Ve arar- om [T freo O at 


=9 [- dt L(r(t), *(2)), (3.3) 


difference: 


where ` represents derivatives with respect to t and the last equality defines the Lagrangian 


Lr): 


pe 


Lir,r) = wP(r) — my] f(r) - FOG (3.4a) 


OF [ap 
Pj= i aah > 0. (3.4b) 


The equations of motion for the particles derived from this action reads 


ay A ae ee TEA 
m m-a A 5) dt Fal f(r) 5) (o 


Since L(r,7) does not depend explicitly on time, its associated Hamiltonian is conserved: 


ay ee 
Or f Ż2 


where 


— uP(r). (3.6) 


This energy represents the difference between the energy of the given configuration and that 
of the infinite unbroken cosmic string. Since the process of nucleation conserves energy, we set 
E = 0. As a result, 7 is determined by the simple condition 


t’ +V(r)=0, (3.7a) 
where 
Vir) = fa ee = i} (3.7b) 


Once can explicitly check that Eqs. (3.7) solve Eq. (3.5). 
The existence of static orbits is equivalent to the requirements 


V(r) = V'(rı)=0. (3.8) 


which readily give 


he s (yP) „ad f(r) = va (VF) PC) 


(3.9) 


rT=riı 


We also want the static orbits to be stable. As usual, stability may be determined by 
computing V”(rı), which turns out to be given by 


v'en =- 2 poy? + EEDE) + alr"). (3.10) 


7 


In particular, if V”(r1) > 0, the static orbit is stable. 

It is interesting to look for solutions defined by spacetimes satisfying the Einstein equations 
sourced by matter that respects the null energy condition. To this end, consider a static 
spherical metric of the form 


dr? 
dst5—=—f(r)d? + — + r'dQ, 3.11 
2D fC ) g(r) 2 ( ) 
where dQ} is the usual round metric on S?, 
dO; = dé? + sin? bdg’. (3.12) 


Consider now a null vector u”. After expressing ut in terms of the other components of the null 
vector u”, contracting u” twice with the matter stress tensor T, will yield 


BTT vuu” = (...)(u")? + (...) [(u’)? + sin? 6(u®)?] . (3.13) 


Requiring the two terms to be non-negative then results in two conditions on the functions f(r) 
and g(r). However, when f(r) = g(r), the first term always vanishes and we are left with 
one condition on f(r). In this case, the non-zero components of the Einstein tensor are 


Gu = po) (—f) (3.14a) 
_ fale pare yt 

Grr = ( 3 ) 7 (3.14b) 
1 G 

Goo = 5r(2f +rf") = Sor (3.14c) 


Using the Einstein equations and noting that the cosmological term vanishes when contracted 
with the null vector gives 


1 
STT tu = (1 -f+ rt") [(u’)? + sin? 0(u®)?] . (3.15) 
The null energy condition then requires 
2(f — 1) 
oe (3.16) 


for all r. 
It would also be reasonable to ask that the stronger weak energy condition holds: T,,,u“u” > 
0 for all timelike observers u”. This gives an additional condition which for f(r) = g(r) takes 
the form: 
3r? 
We will now use what we learned in this section to construct background geometries that 
admit static orbits resulting from the breaking of a cosmic string. 
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4 Semiclassical Production Rate 


Besides finding orbits in a plethora of spacetimes, we are also interested in studying the in- 
stantons that give a finite probability for a test cosmic string to break and produce particles 
attached to the broken ends. To do so we must solve the Euclidean equations of motion and 
find a trajectory r(7) describing expansion from r = 0 to (near) the static orbit r = rı deter- 
mined in the previous section. Here 7 is the Wick rotated time T = it. We use Euclidean time 
translation symmetry to set r = rı at T = 0 and r = 0 at T = 7 < 0. In addition, regularity 
under r + —r demands dr/dr|,=,, = +00 = 7'(0) = 0 (see the blue curves in Fig. |2). 
In the semiclassical approximation, the production rate of such instantons is given by 


rae as (4.1a) 


2 
as--n ff j PBardr +m f Jro T 5 (Z) dr 
(r) aw dr 


= a: uP(r(T)) -mje + a dr. (4.1b) 


where AS is the Euclidean action of the instanton relative to the eternal unbroken cosmic 


string in the background spacetime with S defined to be the analytically continued value of 
the Lorentzian action multiplied by —i, and W is the complement of the Euclidean string 
worldsheet. 

The Euclidean equations of motion for the particles can be readily obtained by varying 
AS with respect to r(r). Just as for Lorentzian signature, these instantons can be reduced to 
studying a single ODE via “energy conservation”, which in this case gives 


7 V | ae (4.22) 


C1 + £P(r)]’ 


where Ci parametrises the “Euclidean energy”. Regularity of T(r) at r = 0 demands C; = 0. 


We thus have 
=4 f r), FY rena (4.2b) 


Inputting the expression for dr(7)/d7 into the on action gives a simple expression for AS 


ere 


Since this expression only depends on 7, through r(7v), we can perform another change of 


(4.3) 


variable and write everything in terms of r instead. This simplifies our final expression to 


(4.4) 


s= f 


This final expression bypasses the need to actually determine the instanton orbits r(7), since 
it depends on f and g only. Note that rı and u/m are determined by Eq. (3.9). 

This is the expression we are going to use in the following sections. For finite r1, it is clear 
that AS is finite. This may at first seem surprising as in the limit where the created particles 
sit on static worldlines r = rı the Euclidean time 7) at which r = 0 must diverge to —oo. 
This is the phenomenon described in the introduction where the Euclidean worldline becomes 
a highly eccentric closed loop that extends much farther in Euclidean time than it does in 
space. In this limit, the action AS receives contributions from an infinite set of Euclidean 
times. However, the contribution from Euclidean times near T = 0 vanish in this limit as any 
static orbit has Euclidean action E f dr, so energy conservation guarantees such contributions 
to precisely cancel analogous contributions to the action of the unbroken cosmic string. The 
net result is that AS receives non-negligible contributions only from Euclidean times near +79 
and thus naturally remains finite as Tọ — —oo. 


5 An AdS Star Example 


In this section we present a simple analytic 3+ 1-dimensional metric satisfying all the conditions 
derived in Sec. 3. In other words, it contains stable static solutions for particles attached to a 
cosmic string which have the same energy as the unbroken string. The metric is: 


dr? 
ds? = —f(r)dt? + —~ + r7dQ2, 5.la 
( ) f(r) 2 ( ) 
2 
1+ - = AA r>R 
M=} Ẹ r (5.1b) 
1+ 55 + Ar’ — Br’, r< R, 


where A and B are positive parameters to be determined, and R is for now an arbitrary 
parameter bigger than the horizon radius for the given M. 

The null energy condition is satisfied everywhere, which can be checked using Eq. (3.16). 
Matching f(r) and f’(r) at R gives 


3M 
= oR ee) 
7M 


The static and energy conservation conditions require the radius of static particles rı to be 


TOKON s 


with rı < R. Thus, R < 3M, i.e. the star has to be fairly compact. 
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If we want the weak energy condition to hold, we also need Too > 0. We find: 
21M(R* — rf) 
4R14 


There is only one term that could be negative, and we can ensure the expression is always 


positive by taking M < i Note that this does not conflict with the previous inequality 


3M > R. In fact, combining these inequalities yields (R/L)? > 7/6, so the star must be larger 
than the AdS radius. Interestingly, note that Too always vanishes as r — R, so that the energy 


87 Too = 


r2 
2 (1 + 5) R” + Mr?(3r* — 7R*)| . (5.5) 


density on the edge of the star is zero. So, in a sense, it is more like a compact spherical cloud 
than a compact spherical star. It is also interesting to note that the weak energy condition fails 
if f(r) = 1 +r?/L? + Art — Br? (for any allowed R and M). This is why we chose an r® term 
in our metric. 
One can explicitly compute the semi-classical production rate in this example. If we define 
x and y by 
2 R? 


=3yM = -— . 
R=3y and M 77a (5.6) 


then the instanton action is given by 


7 1 1 3 2 Ye T 
AS = V2- eve | = arct a4 
Viva al seve |3 E on( 3ye ) 1 
3y? 34/4 + 2ye 2 
ue arctanh V8y4+ 2H + ve (5.7a) 
24/4 + 2ye + y2 4+ ye 


wWIV1—« 1 27x3/2y 
c = ———_ sinh 4 = inh | ——————__ | >. : 
yY avai sin f arcsin Par — A \ (5.7b) 


In particular, the expression is finite as expected. 


where 


6 Vacuum Examples 


We now construct a class of example spacetimes sourced only by a cosmological constant for 
which static, stable orbits can be found, and we compute the associated instanton production 
rate. Global AdS does not satisfy our condition since the Æ = 0 condition would require 
the particles to be created infinitely far apart. But one can view this as being close to satisfying 
the condition, and hope that a small perturbation will allow the particles to be created at a 
finite separation. We will show that this is indeed the case. 

We construct static, axisymmetric solutions for which the boundary metric is not exactly 
the Einstein static universe. We seek bulk geometries for which the boundary metric takes the 
following simple form 

ds? = —dt? + L? [1 + €S,(0)] dn? (6.1) 
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where S,;(@) is a scalar harmonic which, under our symmetry assumptions, is related to the 
Legendre polynomials P; via S;(@) = P;(cos@). To proceed, we will use a mixture of analytical 
and numerical techniques, each of which is explained in a separate subsection. 


6.1 Nonlinear perturbative approach 


In this subsection we will take the so called quasi-spherical gauge, where our bulk spacetime 
metric takes the following form 


2 r? 2 dr? 2402 
ds“ = —F(r,@) {14 T2 dt* +4 H S(r, A)r*dQs , (6.2) 
G 


where F, G and S are to be determined in what follows. For F = G = S = 1 we recover pure 
AdS written in standard global coordinates. We seek solutions for which 
lim F = lim G=1, and limS=1+€S,(0). (6.3) 
Too Too T—> 00O 
We will demand regularity at the poles of the deformed S°, located at 0 = 0,7 and at the 
centre we impose that the geometry is smooth. This in turn implies 


ac 


OF 
G0) = S(0,0), and ar — 


_ as 
2 


Or 


=0. (6.4) 


r=0 r=0 


Since we are interested in this section in a perturbative expansion in €, we shall expand 
all our functions in terms of tensor derived spherical harmonics (see for instance [21]). Given 
that we are working in four spacetime dimensions, there are no tensor harmonics, so we only 
have the scalar derived and vector derived tensor perturbations. The geometries we seek are 
axisymmetric with respect to 0/0¢ and it is a simple exercise to show that no vector harmonics 
on the two-sphere exists within such symmetries. We are thus left with scalar derived tensor 
perturbations. For a generic gauge, and imposing staticity, these are given by 


ogt, = fulr)Se, 659°. = for(r)Se, Ôg, = Hr (r)SGi; + Hr(r)Si;, and 6g‘, = f,(r)WeSe, 

(6.5) 
where i, ïj are indices on the two-sphere, W is the standard Levi-Civita connection on S?, G is 
the round metric on the two-sphere and S,; is a traceless transverse tensor with respect to G, 
given by 


LL +1 
Any linear diffeomorphism, say €, can equally be decomposed in terms of such scalar harmonics 


E = é (r)Sgdr + K(r) (VoSe)d0 P (6.7) 
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Our gauge freedom amounts to choosing ér and K(r) so that Hr(r) = fr(r) = 0. It is a 
simple exercise to show that this gauge choice is possible to achieve, by looking at the relevant 
components of the gauge transformed metric perturbation 


ae 
09a = Öga + V (aS) . (6.8) 


This means, in our gauge, the metric functions can be expanded as 


F(r,0)=1+ Y X ef (r)S.(8), (6.9a) 


i=1 ¢=0 


G(r,0) =14+ X X eg! (r)Se(0), (6.9b) 


S(r,0) =1 + 2 NO €s (r)Se(0). (6.9c) 


We shall present generic results for arbitrary values for even £ to first order in €, and will 
specialise to the higher order case for = 2. To linear order in €, we find 


FP gy gf) 0, for €=Oand/>3, (6.10a) 
and 
gP (r) = AP), ee 
(r= (@+2)T Cay 1 (Jar l-1 T E (6.10¢) 
= Fares yg A er ee T i 
1+ = 
L2 
(1) 
amaa wR f Pa, 2 rna 1 
aa Wr L? (2 + 0 —2) 1+ Bye oD! P40 -IE tae ome 
T 


where F(a, b;c;z) is the Gauss hypergeometric function. The normalisation was chosen so 
that 
lim sP(r)=1. (6.11) 


We are now ready to use Eq. (83.9) to linear order in e. Along 0 = 0,7 and for even £, we 
have 


f(r) = E -+ f(r) (1 + =) + Ole’), (6.12a) 
g(r) = 1 + TRG) (1 -+ =) + O(c"). (6.12b) 
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which gives 
P(r) =r+O(e*). (6.13) 


Once the dust settles, one finds for the critical location of the orbit: 


+ O(e°) (6.14) 


We will want to compare this result to our full nonlinear numerical results of the next 
section. In order to perform such comparison, we need a gauge invariant measure for the 
distance between the two orbits. This is best done by looking at the proper distance between 
the two orbits, given by 


Py =2 


dpi ioe (=) + O€) (6.15) 


vA 1 
i J G(r, 0)4/1+ 5 


The calculation of the on-shell action is very simple to this order in e€, giving 
AS; =T mL + O(e loge) (6.16) 


with the next order term being possible to determine only if F, G and S are known to second 
order in €. 

For particular values of £ we can do substantially better. For instance, we took £ = 2, 
and solved the resulting equations up to second order O(c?) and thus determined the order 
O (eloge) term missing in Eq. (6.16). The intermediate expressions are too cumbersome to be 
reproduced here, but follow mutatis mutandis calculations done elsewhere [24]. 

For the proper distance we now find 


32 384C + 3424 — 457°) L 
Pron = 210g ( Jri = uve: 


2 
FF + O(c), (6.17) 


where C is Catalan’s constant. For the instanton production rate we find 


3 T STE 
ASp=2 = TLM fı + I6 h + 3 + 2 log 63 l +O (eloge) . (6.18) 


Computing AS- to this order in € requires using the following Feynman trick twice: 


ax 
1 + 029242 


A 


1 
arctan(aé) = | dy where a,ĉ€ (6.19) 
0 
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6.2 Exact numerical construction 


We now move beyond perturbation theory and numerically construct exact vacuum solutions 
with deformed boundary metrics. We will use the DeTurck formalism, which was first proposed 
in [25] and reviewed in [26] 27]. Just as the preceding sections, we work in four bulk spacetime 
dimensions, but our methods work equally well in higher dimensiong"] 

The starting point of the DeTurck formalism is to determine the most general ansatz for the 
metric compatible with the desired symmetries. We are interested in static and axisymmetric 
solutions. This means a hypersurface orthogonal Killing field k = 0/0t exists and is everywhere 
timelike as well as a stationary Killing field m = 0/0¢, whose orbits are isomorphic to U(1). 
Using t and ¢ as coordinates leads to the following general line element 


ds? = Gu (r, 0)dt? + Grr, 0)dr? + 2G,6(r, 0)dr dé + Goo(r, 0)d0? + Gyg(r,O)dd?, (6.20) 


which still exhibits full reparameterization invariance in the (r,@) coordinates. The remaining 
gauge freedom is fixed via an appropriate choice of the so called reference metric, which is 
an essential ingredient of the DeTurck formalism. One next imposes the Einstein-DeTurck 
equation: 


3 
Rab + T29% = V (a&b) =0 ) (6.21) 


where €* = g% {.(g) —T,(g)], [%,(g) is the Levi-Civita connection associated with a metric 
g, and g is the so called reference metric. The reference metric is chosen to have the same axis 
location as the metric we seek to find, and satisfies the same Dirichlet boundary conditions as 
g. Solutions of Eq. do not necessarily coincide with solutions of the Einstein equation, 
except if €* = 0. However, it turns out that this is always the case for stationary solutions 
enjoying the so called (t, ¢) symmetry [28], and thus, in particular, for the static solutions we 
wish to determine (for the static case an alternative proof for the equivalence of solutions of 
the Einstein and Einstein-DeTurck equations was presented in [29]). 

The next step is an appropriate choice of reference metric. We motivate it by starting 
in global AdS in standard (t,r,6,@) coordinates and by preforming the following change of 


variables 
yy 2- y’ 
r=L p > cos = rvy2— r?, and t= LT, (6.22) 
=y 
which transforms 
2 r? 2 dr? 2 2 . 2 2 
ds* = — I+ dt“ + 72 +r? (d6? + sin? 8de’) (6.23) 
lcs 
into 
L? 4 dy? A dx? 
eT 2 (Qiang? l=2'fde"|>. 6.24 
ds? = poy {arts SOG tte — EH- aag (6.24) 


4Tn odd spacetime dimensions we would have poor convergence properties because of conformal anomalies. 
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For the reference metric we then take 
p 4 dy? 


re 2 
wga tag 


re PHE) [PS 


where 
H(z,y) = (1-4)? +2- yha), with ~—h(x) = 1 + eS (2v2 = 7) . (6.25b) 


Note that H(x,0) = 1 and H(zx,1) = h(x). Our metric ansatz now reads 


if 4 q(x, 
ds? = f -aenar A (dy + y g(x,y) de)? + 


4qa(z, y) dz? 
2— x? 


y (2—y)H (zy) + (1— 2°)? amg) ae? , (6.26) 
where g(x,y), i € {1,...,5} are functions of x and y to be determined next. We will be 
interested in considering even values of £ only, for which S;(@) is even under the reflection 
0 — m — 0. This translates into parity symmetry under the transformation x — —x, which 
we can use to half our computation domain. So we take (x,y) € [0,1]?, with x = 0 being a 
symmetry hyperplane, y = 0 the regular centre, y = 1 the location of the conformal boundary 
and x = 1 the north pole of the deformed 9°. 

Next we address the issue of boundary conditions. At y = 0, our geometry should locally 
look like flat space, because this corresponds to a point. This in turn enforces 


Oqi(@, y) 


qalx, 0) T qs(2, 0) B qo(x, 0) and By = 0. (6.27) 
Y ho 
At x = 0 reflection symmetry imposes 
dq (x,y) Oqo(x, y) Oqa(x, y) ðqs(x, y) 
= = = =0 d 0,y)=0. 
Ox x=0 Ox x=0 Ox x=0 Ox x=0 = o 4) 
(6.28) 
At x = 1, the north pole of the deformed two sphere, smoothness imposes 
On(z,y)| _ aleu) _ ðalzy)| _ gley) o 
Ox x=1 Ox x=1 Ox x=1 Ox r= 
q3(1, y) — 0, and qa(1, y) = qs(1, y) : (6.29) 


Finally, at the conformal boundary, we take q(x,1) = qo(z,1) = qi(z,1) = gs(z,1) = 1 
and q3(x,1) = 0, which means that the line element (6.26) is chosen to approach the reference 
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metric (6.25ah. All we need to show is that the reference metric has the correct boundary 
metric. This is easily seen if we use the coordinate transformation (6.22). 

To solve the Einstein-DeTruck equation numerically, we discretise the equations using a 
Chebyshev grid on Gauss-Lobatto nodes and solve the resulting equations using a standard 
Newton-Raphson routine. Such methods have been reviewed extensively in [27]. 


6.3 Results and comparisons 


Having constructed the new vacuum solutions, we now ask whether there are stable, static 
solutions for particles attached to cosmic strings with the same energy as the unbroken string. 
We first note that the deformation parameter, €, is necessarily bounded above, since the bound- 
ary metric would otherwise have a curvature singularity (this is best illustrated by computing 
the Ricci scalar of the boundary metric) at the location of the zero of 1 + €S;(@). In Fig. 
we plot €max(Z) as a function of £ up to 4 = 20. It is easy to show that, asymptotically, 
Ga = + O(1/0), where J,(x) is a J Bessel function of order k, and jg; is the / zero 
of the J Bessel function of order k. 


2.4 


2.3- 


Emax (4) 


2.11 


2.0 


Figure 3: €max as a function of £: the blue disks are the exact numerical data and the red dashed 
curve the leading approximation at large £. 


We begin by studying the case with £ = 2, for which €,,,, = 2, and briefly mention at the 
end what happens for different even values of l > 2. First, we would like to understand whether 
stable orbits exist for arbitrarily small values of €, as predicted by our perturbative non-linear 
results. This indeed turns out to be the case. In Fig. [4} we compute the proper distance between 
the two orbits as a function of e for l = 2. The blue disks correspond to our exact numerical 
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results and the solid red line to the analytic expression (6.17): the agreement between the two 
when e < 1 is reassuring. As we increase the deformation on the boundary, the separation 


0.0 0.5 1.0 1.5 2.0 


Figure 4: Proper distance between the two stable orbits computed for £ = 2: the blue disks 
correspond to our exact numerical data and the solid red line to the perturbative expression 
(6.17). 


between the static particles decreases, but always stays larger than the AdS radius. 

Next we compute the instanton action by computing using our numerical solutions. 
The result of this calculation is seen in Fig. |5| where our numerical data is represented by 
the blue disks and the solid red line corresponds to our perturbative result (6.18). Again, the 
agreement between these two expressions at small values of e is indicative of the correctness 
of both methods. As expected, the action decreases as the separation between the particles 
decreases. Note that the nonlinear perturbative solution and the exact numerical solution are 
obtained in different gauges. 

We have extended our calculations to other even values of £, and find that increasing £ 
increases P; for large € and decreases P; for small enough e. We can observe this effect on the 
left panel of Fig. [6] where we plot the proper distance between orbits for several values of £, 
each of which are labelled in the figure caption. The behaviour at small € is corroborated by 
the nonlinear perturbative prediction (6.15). Nevertheless, for each value of £, we observe that 
increasing € does decrease Py. 

The behaviour of the instanton action for various £ is shown on the right panel of Fig. [6] 
The action increases with @, and unlike the £ = 2 case, for sufficiently large Z the action grows 
with e. 
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3.0 
2.8 
y 
ASi=2 
mL z6 
2.4 
2.2 
0.0 0.5 “1.0 15 2.0 


Figure 5: Instanton action computed for the £ = 2 profile: the blue disks correspond to our 
exact numerical data and the solid red line to the perturbative expression (6.18). 


7 Discussion 


We have presented examples of spacetimes in which a (test) cosmic string can break and produce 
two particles with arbitrarily small acceleration. There is a long history of showing that small 
black holes follow geodesics (see [30] for a fairly recent discussion and references to earlier work). 
Similar arguments should show that particles attached to cosmic strings can be replaced by 
small black holes with no change in their dynamics. (The Euclidean solutions may have conical 
singularities on the horizon but they still describe pair creation .) 

Previous examples of black hole pair creation show that the black holes are created with 
their horizons identified at the moment of creation [10] and therefore that the instanton creates 
a wormhole. In the limiting case of zero acceleration and for instantons free of conical singulari- 
ties, the path integral computing the wormhole’s quantum state is precisely of the form studied 
in [6], so with appropriate boundary conditions the back-reaction of quantum fields will render 
the wormhole traversable in the semiclassical approximation. (We again remind the reader 
that some large parameter is generally required to render quantum fluctuations small [5].) As 
argued in the introduction, the same will be true with sufficiently small non-zero acceleration; 
see in particular for the condition relevant to the traversable wormholes of [5]. 

The net result is the creation of a traversable wormhole by a nonperturbative process in 
quantum gravity. We have seen that while this is not possible using only light cosmic strings 
in pure AdS, it is possible in nearby vacuum solutions. It is interesting to note that these 
solutions all have negative total energy [BI]. This leads to the following intuitive picture: the 
attractive force of pure AdS is not quite sufficient to create the desired trajectories, but when 
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3.21 


3.0| £ 

| -2 

AS agl = 4 
mL + 

+6 
2.6} 

t a 8 

24! æ 10 


ĝo o5 “To 15 20 260 05 «0 15 20 


Figure 6: Proper distance between orbits (left panel) and instanton production rate (right 
panel) computed for ¢ = 2,4,6,8,10 plotted as a function of e, with each curve labelling a 
different value of £. 


we deform the boundary metric we lower the energy, effectively producing a deeper potential 
well in which it is possible. 
When we replace the test particle with a small black hole, we only change the Euclidean 


solution in a small tubular neighborhood of the particle trajectory which has topology S1 x R3. 
The S! represents the Euclidean time along the particle worldline. To have smooth instantons, 
this period must match the inverse temperature of the black hole. Since the particles in our 
instantons travel a long Euclidean proper distance before their orbits close, we can find such 
smooth black hole instantons by pasting in nearly-extremal Reissner-Nordstrom or Kerr black 
holes. If the right sort of fields propagate in the wormhole throat, this should result in a 
wormhole that is traversable for long periods of time. The perturbative argument is given in [6], 
and it is also expected from continuity with the non-perturbative constructions of [4] [5] which 
are related to thermofield double states at extremely small temperatures. If the acceleration is 
small enough, we may expect the instanton to produce wormholes that remain open forever. 
However, one should recall that instantons with conical singularities can sometimes yield 
higher pair-creation rates than smooth ones (see e.g. [19]). This seems likely to occur in our 
context, as black holes farther from extremality will have greater entropy, so reducing the charge 
and angular momentum at fixed mass should increase the rate of pair production} If one seeks 


>In general, smooth instantons will extremize the action subject to varying the black hole area while holding 
fixed other parameters and boundary conditions; see e.g. the recent discussion in [32]. But in [19] this extremum 
has a negative mode and is not a minimum. In addition, in our case we remove the singularity by adjusting not 
only the black hole area, but also both the charge q and associated potentials at infinity while holding fixed the 
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a traversable wormhole that remains open for long times, it may thus not suffice to simply wait 
for the first wormhole to appear on our static trajectory. One must also wait some additional 
time for one that is nearly extremal. On the other hand, the fact that energy conservation fixes 
the mass m of any particle created by our instanton near the static orbit leads to a positive 
feature of our process: one may choose parameters (including the string tension u) to make 
this m large, perhaps even larger than our current approximations justify, and then have the 
comfort of knowing that only large wormholes will be created at that location] 

It is not yet clear how large is the class of spacetimes which allow the creation of a traversable 
wormholes. In the vacuum examples discussed in the previous section, the particles are always 
created with a separation larger than the AdS radius L. We have also tried many other 
deformations of the boundary metric to see if we could lower the distance between the two 
particles below L, but were unable to do so. This suggests that the AdS attraction is playing a 
key role and similar phenomena may not be possible in vacuum asymptotically flat spacetimes. 

In addition to the examples discussed here, we have also found further examples by deform- 
ing the AdS soliton [83]. However, in that case there was a minimum deformation required to 
find stable static orbits. Also in that case, we were never able to decrease the distance between 
the two orbits to be below the AdS length scale L. 

We also tried a number of asymptotically AdS spacetimes for which it does not seem to be 
possible to nucleate a traversable wormhole. We have constructed (perturbatively) a number 
of boson stars in AdS for a wide range of mass parameter, with alternative and standard 
quantization, and we were not able to find stable orbits. We also constructed smooth geometries 
in global AdS with scalars satisfying double trace boundary conditions, and for those, also no 
orbits were found. 
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location of the trajectory. Extremizing the action under such variations need not remove conical singularities. 

Replacing our cosmic string with a magnetic flux tube should lead to an interesting variation. At first order 
in the magnetic charge q, the mass of particles created by decay of the magnetic field at a given separation is 
proportional to q. So watching a fixed static orbit selects a given charge-to-mass ratio q/m rather than the fixed 
mass m obtained in our case. One should thus be able to engineer an instanton that creates only nearly-extreme 
black hole pairs regardless of whether conical singularities are allowed at the horizon. But then the dominant 
wormholes created on the static orbit are likely to be small and large ones will be more rare. 
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